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Abstract 



The cooling of a compact star depends very sensitively on the state of dense matter 
at supranuclear densities, which essentially controls the neutrino emission, as well 
^ . as on the structure of the stellar outer layers which control the photon emission. 

\ Open issues concern the hyperon population, the presence of meson condensates, 

' superfluidity and superconductivity, and the transition of confined hadronic matter 

■ to quark matter. This paper describes these issues and presents cooling calculations 

, based on a broad collection of equations of state for neutron star matter and strange 

^ ' matter. These results are tested against the body of observed cooling data. 

Qh| Key words: Nuclear matter, Quark Matter, Equation of state. Neutron stars, 

Q ■ Cooling 

PACS: 97.10.Cv, 97.60.Jd, 26.60+c, 12.38.Mh 



1 Introduction 



Cooling simulations, confronted with soft X-ray, extreme UV, UV and optical 
observations of the thermal photon flux emitted from the surfaces of neutron 
stars, provide most valuable information about the physical processes operat- 
ing in the interior of these objects. The predominant cooling mechanism of hot 
(temperatures T ^ 10^° K) newly formed neutron stars is neutrino emission, 
with an initial cooling time scale of seconds. Neutrino cooling still dominates 
for at least the first thousand years, and typically for much longer in slow 
(standard) cooling scenarios. Photon emission eventually overtakes neutrinos 
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when the internal temperature has sufficiently dropped. Being sensitive to the 
adopted nuclear equation of state (EOS), the stellar mass, the assumed mag- 
netic field strength, superfiuidity, meson condensates, and the possible pres- 
ence of color-superconducting quark matter, theoretical cooling calculations 
serve as a principal window on the properties of super-dense hadronic mat- 
ter and neutron star structure. The thermal evolution of neutron stars also 
yields information about such temperature-sensitive properties as transport 
coefficients, crust solidification, and internal pulsar heating mechanisms. 

We will present here an overview of the current status of neutron star cool- 
ing calculations tested against the steadily growing body of observed cooling 
data on neutron stars. The reader can also find a complementary approach in 
Ref. [1]. Space limitation forbid us to dwell into the discussion of observational 
data and we refer the reader to Ref. [2] for a recent compilation and to two 
recent observational reviews [3,4] to get a flavor of the difficulties involved in 
the analysis and interpretation of the data. 

The paper is organized as follows. In § 2 we introduce the basic equations and 
the physics input that governs neutron star cooling. Some simple analytical 
solutions to the cooling equations are presented in § 3. The minimal cooling 
paradigm, which assumes that no enhanced neutrino emission is allowed in 
neutron stars, is presented in § 4. In § 5 we discuss enhanced neutron star 
cooling via the direct Urea process, meson condensates, and quarks. The cool- 
ing behavior of compact stars made of absolutely stable strange quark matter 
is explored in § 6. The impact of magnetic fields on cooling and heating mech- 
anisms are discussed in § 7 and § 8, respectively while § 9 considers cooling 
neutron stars in transiently accreting binary systems. Conclusions are offered 
in 5 10. 



2 Basic Equations and Physics Input 

The basic features of the cooling of a neutron star are easily grasped by simply 
considering the energy conservation equation for the star. In its Newtonian 
formulation this equation reads 



where i?th is the thermal energy content of the star, T its internal temperature, 
and Cv its total specific heat. The energy sinks are the total neutrino luminos- 
ity Lj,, described in § 2.2, and the surface photon luminosity L^, discussed in 
§2.5. The source term H includes all possible "heating mechanisms" which. 





dt 
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for instance, convert magnetic or rotational energy into heat as summarized 
in § 8. Some simple analytical solutions to Eq. (1) will be presented in § 3. 

The dominant contributions to Cy come from the core, constituting more 
than 90% of the total volume, whose constituents are quantum liquids of lep- 

tons, baryons, mesons, and possibly deconfincd quarks at the highest densities. 
When baryons, and quarks, become paired, as briefly described in § 2.3, their 
contribution to Cy is strongly suppressed at temperatures T <C Tc (T^ being 
the corresponding critical temperature). The crustal contribution is in prin- 
ciple dominated by the free neutrons of the inner crust but, since these are 
certainly extensively paired, practically only the nuclear lattice and electrons 
contribute. Extensive baryon, and quark, pairing can thus signiflcantly reduce 
Cy, but by at most a factor of order ten since the leptons do not pair. 

All results presented in this chapter were obtained with numerical codes which 

exactly solve both the energy conservation and the heat transport equations 
in their General Relativistic forms (see, e.g., [2,5]). Since the density and 
chemical composition of the star, after the proto-neutron star phase, do not 
change with time, the Tolman-Oppenheimer-Volkoff equation of hydrostatic 
equilibrium is solved initially, for a given EOS, and only the thermal equations 
are evolved with time. The outermost, low density, layers of the star do see 
their structure evolve during the cooling and they are hence treated separately 
as an envelope which is the outer boundary condition and is described in § 2.5. 



2.1 The Equation of State (EOS) 



The cross section of a neutron star can be split roughly into three, possibly 
four, distinct regimes. The first regime is the star's outer crust, which consists 
of a lattice of atomic nuclei and a Fermi liquid of relativistic, degenerate elec- 
trons. The second regime, known as the inner crust and where free neutrons 
appear, extends from neutron drip density, ~ 4 — 8 x 10^^ g cm~^, to a transi- 
tion density of about 2 x 10^"^ g cm~^. Beyond that density one enters the star's 
third regime, its core matter where all atomic nuclei have dissolved into their 
constituents, protons and neutrons. Furthermore, because of the high Fermi 
pressure, the inner core may be expected to contain baryon resonances, boson 
condensates, hyperons, or/and a mixture of deconflned up, down and strange 
quarks [6,7]. The EOS of the outer and inner crust is rather well known. This 
is very different for the EOS of the star's core which is only very poorly un- 
derstood. Models derived for it fall into non-relativistic EOS and relativistic, 
field theoretical ones. 

The most frequently used non-relativistic models for the EOS are based on 
the hole-line expansion (Brueckner theory), the coupled cluster method, self- 



3 




e (MeV/fm') 

Fig. 1. Models for the EOS of high-density neutron star matter computed for dif- 
ferent compositions and many-body techniques. (p,n denote protons and neutrons; 
H, K~, Q stand for hyperons, K~ condensate, and quarks, respectively.) 

consistent Green functions, the variational approach (Monte Carlo techniques), 
the semi-classical Thomas-Fermi method, and the density functional approach 
based on Skyrme effective interactions (for an overview of these methods 
and additional references, see for instance Refs. [8,9,10]). The forces between 
hadrons are described in terms of phenomenological nucleon-nuclcon inter- 
actions, possibly supplemented with three-nucleon interactions to achieve a 
better agreement with the binding energy of nuclear matter at the empirical 
saturation density. Figure 1 shows two sample models for the EOS of neutron 
star matter (protons and neutrons only) which arc obtained from variational 
calculations based on the Urbana Vu two-nuclcon interaction supplemented 
with the three-body interactions UVII (left curve) and TNI (right curve) [11]. 

Relativistic, field-theoretical approaches to the EOS are based on model La- 
grangians which generally describe baryons as Dirac particles interacting via 
the exchange of mesons. The most important mesons are the a and the uj, 
which are responsible for nuclear binding, while the p meson is required to ob- 
tain the correct value for the empirical symmetry energy. Nonlinear a terms 
need to be included at the mean-field (relativistic Hartree, RH) level in order 
to obtain the empirical incompressibility of nuclear matter [6]. Such terms arc 
not necessary if the field equations are solved for the relativistic Brueckncr- 
Hartree-Fock (RBHF) approximation [7] . In contrast to RH, where the param- 
eters of the theory are adjusted to the properties of infinite nuclear matter, 
the RBHF method makes use of one-boson-exchange potentials whose param- 
eters are adjusted to the properties of the deuteron and to the nucleon-nucleon 
scattering data. Figure 1 shows several sample EOSs computed for RH and 
RBHF, assuming different particle compositions of neutron star matter [7]. 
Finally we mention that in recent years a new class of effective field theories 
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Fig. 2. Neutron star mass versus radius for different EOSs. The horizontal lines 

refer to the masses of Vela X-1 (1.88 ± 0.13 Mq) [15], Cyg X-2 (1.78 ± 0.23 M©) 
[16], and PSR 1913+16 (1.4408 ± 0.0003 M©) [17] (in all three cases quoted mass 
uncertainties are la). 

was developed which treat the meson-micleon couplings as density dependent. 
These field theories provide a very good description of the properties of nuclear 
matter, atomic nuclei as well as neutron stars [12,13,14]. 

In Figure 2 we show the mass-radius relationships of neutron stars for different 
EOSs: it illustrates the well-known fact that this relationship is very sensitive 
to the underlying model for the EOS. 

Phase transitions to boson condensates, hyperonic matter, or quark matter, 
soften the EOS, hence reducing both the stars gravitational mass and radius 
and resulting in a smaller maximum mass sustainable by the EOS. A generi- 
cally different mass-radius relationship is obtained for stars made of absolutely 
stable strange quark matter (strange stars). For them M oc R^, which is only 
significantly modified if M is close to the mass peak. As pointed out in [18], 
strange stars can carry a solid nuclear crust whose density at its base is strictly 
limited by neutron drip. This is made possible by the displacement of elec- 
trons at the surface of strange matter, which leads to a strong electric dipole 
layer there. The associated electric field is so strong that it holds open a gap 
between the nuclear crust and quark matter, preventing the conversion of the 
crust into the hypothetical lower-lying ground state of strange matter. Obvi- 
ously, free neutrons, being electrically charge neutral, cannot exist in the crust 
because they do not feel the Coulomb barrier and thus would gravitate toward 
the strange quark matter core where they are converted, by hypothesis, into 
strange matter. Consequently, the density at the base of the crust (inner crust 
density) will always be smaller than neutron drip. The situation is graphically 
illustrated in Figure 3 which shows the EOS of a strange star [7,19] carrying 
a nuclear crust. Since the crust is bound to the quark matter core by gravity 
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Fig. 3. EOS of strange quark matter surrounded by nuclear matter. The maxi- 
mal possible nuclear matter density is determined by neutron drip which occurs at 
Ccrust = 0.24 MeV fm~^ (4.3 x lO""^^ g cm~^). Any nuclear density that is smaller 
than that is possible. As an example, we show the EOS for a sample density of 
ecrust = 10-^ MeV fm-^ (10^ g cm-^). 

rather than confinement, the mass-radius relationship of a strange star with a 
nuclear crust is qualitatively similar to the one of purely gravitationally bound 
neutron stars, as shown in Figure 2. 

2.2 Neutrino Processes 

As already mentioned in the introduction, neutron stars are born with temper- 
atures in excess of lO^'^ K. The dominating cooling mechanism of such objects 
for the first several thousand years after birth is neutrino emission from the 
interior. After that, cooling via photon emission from the star's surface takes 
over. Tables 1 and 2 summarize the dominant neutrino emitting processes to- 
gether with their efficiency for neutron star cooling and we now briefly describe 
the most important ones. The reader is referred to Refs. [7,20,21,22] for more 
details. We will nevertheless comment in § 2.6 on minor neutrino emission 
processes which may become important in cases where all the dominant ones, 
as discussed in § 2.3, are suppressed by pairing. 

Direct Urea Processes. The beta decay and electron capture processes 
among nucleons, n ^ p-\- e~ -\- D^, and p + e~ n + i/g, also known as nucleon 
direct Urea process (or cycle), are only possible in neutron stars if the proton 
fraction exceeds a critical threshold [23]. Otherwise energy and momentum 
can not be conserved simultaneously for these reactions. For a neutron star 
made up of only neutrons, protons and electrons, the critical proton fraction 
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Table 1 

Dominant neutrino emitting processes in neutron star cores, in absence of hyperons" 
and quarks^. 



Name 



Process*^ 



Emissivity'' 



(erg cm s 



3 o-i^ 



Modified Urea cycle 
(neutron branch) 

Modified Urea cycle 
(proton branch) 

Bremsstrahlung 

Cooper pair 
formations 

Direct Urea cycle 

7r~ condensate 
K~ condensate 



?T- + p + e^ — >n + n + fe 

p + p + e" p + n + Ue 

n+n^n+n+v+v 
n+p^n+p+v+u 
p + p P + P + u + u 

n + n ^ [nn] + u + v 
p + p ^ \pp] + v + 1/ 

n ^ p + e~ + Ue 
p + e~ n + Ve 

n+ < Tr~ >— > n + e~ + i^e 

n+ < K~ >^ n + e~ + i>e 



2x10^1 R 



10^1 R 



10^^ RT^ 



5x10^^ R Tj 
5x10^9 R 



~ 10'^^ RT^ 
~ 10^5 R T| 



Slow 
Slow 

Slow 

Slow 

Fast 

Fast 
Fast 



" In the presence of hyperons, most processes listed here have corresponding ones 
with hyperons replacing nucleons. 
See Table 2 for quark processes. 

Where are present, all processes involving e~ have a corresponding one with 
li~ , u^, and replacing e~ , Ve, and u^, respectively. 

Quoted emissivities are only indicative: each process has its specific dependences 
on medium and particle densities, effective masses, plus medium effect corrections, 
for which we refer the reader to Rcfs. [7,20,21,22]; the many R factors are the 
respective temperature dependent control functions which take into account the 
effects of pairing as discussed in § 2.3 and detailed in Ref. [21]. 



is around 11%. This follows readily from k^^ = kp^ + kp^ combined with 
the condition of electric charge neutrality of neutron star matter. The tri- 
angle inequality then requires for the magnitudes of the particle momenta 
< + kp^, and charge neutrality constrains the particle Fermi momenta 
according to kp = kp^. Substituting kp = kp^ into the triangle inequality 
leads to kp^ < 2kp^ so that n„ < Sn^ for the number densities of neutrons and 
protons. Expressed as a fraction of the system's total baryon number density, 
n = np + Un, one thus arrives at Up/n > 1/9 ~ 0.11 as quoted above. Medium 
effects and interactions among the particles modify this value only slightly 
but the presence of muons raise it to about 0.15. Hyperons, which may exist 
in neutron star matter rather abundantly, produce neutrinos via the direct 
Urea process E~ — > A + e~ + with A + e" — > E~ + i/e and similar ones 
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involving hyperons and nucleons simultaneously [24,25]. The cooling driven 
by the nucleon direct Urea, however, dominates over the energy loss produced 
by the direct Urea process among hyperons [26]. 

Meson Condensate Urea Processes. The pion or kaon meson fields may 
develop condensates in dense neutron star matter. These condensates would 
have two important effects on neutron stars. Firstly, they would soften the 
EOS above the critical density for onset of condensation, which reduces the 
maximal possible neutron star mass. Secondly, since the condensate, <7r~> or 
<K~>, can absorb as little or as much momentum as required by the scatter- 
ing processes n+ <7r~>— > n+e~ + v'e or n+ <K^>^ n + e~ + Ue, the resulting 
neutrino emissivitics of meson-condensed matter [27,28,29], even though not 
as high as the ones of the direct Urea processes (see Table 1), still lead to 
fast cooling [30,31,32]. Since the K~-condensate process involves strangeness 
change it is less efficient that the 7r~-condensate process, roughly by a factor 
sin^ 9c — 1/20 {9c being the Cabibbo angle). However, medium effects can 
reduce the 7r~-condensate process by about one order of magnitude and make 
it comparable to the i^~-condensate one [33]. Estimates predict the onset of 
charged pion condensation at a density n^j. ~ 2no (no = 0.16 fm~^ being the 
empirical nuclear matter density). However, this value is very sensitive to the 
strength of the effective nucleon particle-hole repulsion in the isospin-1, spin- 
1 channel, which tends to suppress the vr-condensation mechanism and may 
push n^j, to much higher vahics. Similarly, depending on the nuclear model, 
the threshold density for the onset of kaon condensation, n^, is at least of the 
order of 4no [34] . 

Modified Urea Processes. In absence of hyperons or meson condensates, 
or in case the proton fraction is below threshold, none of the above described 
Urea processes is possible. In this case, the dominant neutrino emission process 
is a second order process, variant of the direct Urea process, called modified 
Urea process [35,36], in which a bystander neutron or proton participates to 
allow momentum conservation (see Table 1). Since this modified Urea process 
involves 5 degenerate fermions, instead of three for the direct Urea and meson 
Urea processes, its efficiency is reduced, simply by phase space limitation, by 
a factor of order {T/Tpf. This reduction, for Tp ~ 100 MeV and T = 0.1 MeV 
~ 10^ K, amounts to about 6 order of magnitude (!) and an overall temperature 
dependence instead of T^. It is certainly the dominant process for not to 
high densities in absence of pairing, and is the essence of the old "Standard 
Cooling Scenario" . However, in presence of pairing, neutrino emission by the 
constant formation of Cooper pairs, see § 2.3, most probably dominates over 
the modified Urea process. 

Since the modified Urea process involves a strong interaction for the momen- 
tum exchange between the neutrino emitting nucleon and the bystander one, 
it is prone to medium corrections which seem to result in a reduction of emis- 
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Table 2 

Dominant neutrino emitting processes in deconfined quark matter. 



Name 



Process" 



Emissivity 



(erg cm ^ s 



3 o-l^ 



Efficiency 



Direct Urea cycle 
{ud branch) 

Direct Urea cycle 
{us branch) 

Modified Urea cycle 
(ud branch) 

Modified Urea cycle 
(us branch) 

Bremsstrahlungs 

Cooper pair 
formations 



u + e ^ d + 
d ^ u + e~ + i>e 

li + e~ — J- s + fe 
s ^ u + e~ + Ve 

Q + U + e~ Q + d + Ue 
Q + d^Q + u + e~ + i?e 

Q + u + e- ^Q + s + i^e 
Q + s ^ Q + u + e^Ue 

Qi + Q2^ Q1 + Q2 + V + v 

■u + tt — > [uu] + V + 9 
d+d^ [dd] + U + U 
s + s ^ \ss] + u + u 



10^6 R r| Fast 



10^5 R r| Fast 



1021 R r| Slow 



lO^o R Tl Slow 



~ 10^9 R Tl Slow 

~ 2.5x1020 iiTj 

~ 1.5x10^1 i?rj Slow 

~ 1.5x10^1 i?rj 



" Muons are never present in quark matter, but in case e"*" are present instead of 
e~, processes involving e~ are replaced by similar ones involving e"*". 

Quoted emissivities are only indicative: each process has its specific dependences on 
medium and particle densities, strange quark mass m^, and color coupling constant 
see Ref. [7,40,41] for more details; the many R factors are the respective control 
functions which take into account the effects of pairing as discussed in § 2.3. 

sivity (see. e.g., [37,38]). However, softening of the pion mode, which even- 
tually leads to 7r~-condensation, do results in a very strong enhancement of 
the emissivity when the density approaches n1^, and gives a smooth transition 
from the modified Urea process toward the 7r~-condensate process through a 
medium-modified- Urea process ("MMU" process [22,39]). 

The Quark Urea Processes. The neutrino emission processes in non-super- 
conducting quark matter [40] can be divided into slow and fast ones, in com- 
plete analogy to the nucleon processes discussed above, and the dominant pro- 
cesses are listed in Table 2. The fast quark direct Urea processes d — * u+e'+Ue 
and s u + e~ + Ue are only possible if the Fermi momenta of quarks and 
electrons obey the triangle inequalities associated with these reactions, which 
are kp^^ < kp^ + kp^ and kp^ < kp^-\-kp^. If the electron Fermi momentum, kp^, 
in quark matter is too small for the quark triangle inequalities to be fulfilled, 
a bystander quark is needed to ensure energy and momentum conservation in 
the scattering process. This latter case is referred to as quark modified Urea 
process. The emissivities associated with the quark modified Urea processes 
are considerably smaller than those of the direct Urea processes because of 
the different phase spaces associated with two-quark scattering and quark 
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decay. In the extreme case when the electron fraction vanishes entirely in 
quark matter, both the quark direct and the quark modified Urea processes 
become unimportant. The neutrino emission from the star is then dominated 
by bremsstrahlung processes, Qi + Q2 ^ Qi + Q2 + ^ + ^, where Qi and Q2 
denote any pair of quark flavors. In this case, stellar cooling proceeds rather 
slowly. 

Notice that non- Fermi liquid corrections lead to an enhancement of the emis- 
sivities compared to the results of [40] and Table 2, which for the direct Urea 
is of order {\og{mc^ / ksT))^ , with mc^ ~ 400 MeV [42]. However this effect 
has not yet been incorporated in numerical calculations and is balanced by a 
similar increase in the specific heat so that the net result of non-Fermi liquid 
effect is not expected to be very large [42]. 

2.3 Pairing 

Pairing will unavoidably occur in a degenerate Fermi system in case there is 
any attractive interaction between the particles. In case of the baryons, and 
quarks, in the neutron star interior there are many candidates for channels 
of attractive interactions, and the real question is rather what is the critical 
temperature Tc at which pairing occurs? Calculation of Tc are notoriously 
difficult and results often highly uncertain. We refer the reader to [43,44] 
for detailed reviews. With respect to leptons, there is no obvious attractive 
interaction which could lead to pairing with a Tc of significant value. 

In the case of nucleons, at low Fermi momenta, pairing is predicted to occur 
in the ^Sq angular momentum state while at larger momenta neutrons are 
possibly paired in the — ^F2 state (the mixing being due to the tensor 
interaction). In the case of the neutron ^Sq pairing, which occurs at densities 
corresponding to the crust and, possibly, the outermost part of the core, much 
efforts have been invested in its study and calculations are converging with 
time when more and more sophisticated many-body models are used. In the 
case of the proton ^Sq pairing the situation is more delicate since it occurs 
at densities (in the outer core) where protons are mixed, to a small amount, 
with neutrons. Predictions for Tc span a much wider range than in the case of 
the neutron ^5*0 gap. Whether or not neutrons pair in the ^P2 — channel 
is still uncertain, since even the best available models for the nucleon-nucleon 
interaction fail to reproduce the measured ^P2 phase shift in vacuum [45]. 
Moreover, the results of [46] , which consider polarization contributions to the 
effective interaction, indicate that this gap may be vanishingly small. A set of 
representatives predictions for the nucleon gaps are shown in Figure 4. 

The enormous impact of pairing on the cooling comes directly from the ap- 
pearance of the energy gap A at the Fermi surface which leads to a suppression 
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Fig. 4. Predictions for the nucleon pairing T^. Left panel: neutron pairing; all 
calculations include medium polarization in some measure except for the dotted 
curve which illustrates the strong reduction effect of polarization on this gap. Right 
panel: same for protons; no polarization effects are taken into account except for 
the continuous curve; the shaded area shows the range of expected T^s when po- 
larization is taken into account, according to the estimates of [49]. Central panel: 
neutron pairing, three typical results illustrating the possible range according 
to [45]. See [2] for details and references. 

of all processes involving single particle excitations of the paired species. When 
T <C Tc the suppression is of the order of e~^l^^'^ and hence dramatic. The 
suppression depends on the temperature dependence of A and the details of 
the phase space involved in each specific process. In numerical calculations it 
is introduced as a control function. For the specific heat one has 

Cv(T) cf\T) = Rc{T/T,) X c™'(T) , (2) 



and the control functions have been calculated for both ^5*0 and pairing 
in [47]. For neutrino processes there is a long family of control functions for 
all processes which must also consider which of the participating baryons are 
paired. As for Cy one uses 

e,(T) er^<^(T) = K{T/%) x e™'(T) , (3) 



and the RyS for many processes can be found in [48]. We plot in Figure 5 
several examples of control functions. 

It is important to notice that the gap A is actually a function of the particle 
momentum, A = A(k). For ^5*0 pairing it is isotropic, i.e., A = A(A;), but for 
— ^-^2 pairing the angular dependence of A(k) is complicated. In this latter 
case many phases, with distinct angular dependences of A(k) are possible 
([50] found there are at least 13 of them) and for several of them A(k) has 
nodes at some points of along some lines on the Fermi surface. The control 
functions plotted in Figure 5 assume nodeless gaps, but in cases of ID nodes 
R ~ (T/Tc)^, while for 2D nodes R ~ T /Tc, aXT <^Tc instead of a Boltzmann- 
like suppression. 
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TAc TA, TAc 

Fig. 5. Control functions for Cy (left panel), Eq. (2), 6,^ of the modified Urea process 
(central panel), Eq. (3), and the PBF process, Eq. (5) (right panel). 

Cooper Pair Breaking and Formation (PBF) Processes. Besides the 
above described, and well known, suppressing effects on the specific heat and 
neutrino emissivities, the onset of pairing also opens new channels for neu- 
trino emission. The supcrfluid or superconducting condensate is in thermal 
equilibrium with the single particle ("broken pairs") excitations and there is 
continuous formation and breaking of Cooper pairs, which are very intense at 
temperatures slightly below The formation of a Cooper pair liberates an 
energy which can be taken away by a u-l/ pair [51,52] 

X + X — ^ [XX] + U + V (4) 



where [XX] denotes a Cooper pair of particles X {X stands for neutrons, 
protons, hyperons, quarks, etc.). As an example, the emissivity for neutron 
pairing is [48] ^ 

C"-=8.6xlO='(^)'"(^)x^,,(r/Td(j5^)'. (6) 



where the control function R is plotted in the right panel of Figure 5: the 
process turns on at T = Tc, with an increasing efficiency when T decreases, 
since the energy of the emitted neutrinos is determined by the gap's size 
which grows with decreasing temperature just below T^, and is eventually 
exponentially suppressed when T -C Tc as pair breaking is frozen because 
ksT <C A. This process can be seen as a bremsstrahlung with a very strong 
correlation in the final state and referring to Tables 1 and 2, one sees that 
it is much more efficient that the simple bremsstrahlung one and it can even 
dominate over the standard modified Urea process. This process, which is 
as standard as the modified Urea process, is an essential ingredient of the 
Minimal Cooling Paradigm described in § 4. Analogous processes occur for all 
cases of pairing: neutron, proton, hyperons, and quarks [41]. 



^ Our control function R differs from the R in [48] in that it is normalized to have 
a maximum value of one. 
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2-4 Color- Superconductivity of Quark Matter 



Already several decades ago it had been suggested that the attractive force 
among quarks may cause them to form Cooper pairs [53,54]. Originally the 
gap was estimated to be around A = 0.4 MeV [54]. Recently, however, it was 
discovered that the condensation patterns of quark matter arc much more 
complex than originally thought [55,56]. This has its origin in the fact that 
quarks come in different colors, different flavors, and different masses. More- 
over, bulk matter is neutral with respect to both electric and color charge, and 
is in chemical equilibration under the weak interaction processes that turn one 
quark flavor into another. Depending on density and temperature, quarks may 
thus condense in one of the following pairing schemes. At asymptotic densi- 
ties the ground state of (3-flavor) QCD with a vanishing strange quark mass 
is the color-flavor locked (CFL) phase in which all three quark flavors par- 
ticipate symmetrically. This phase has been shown to be electrically neutral 
without any need for electrons for a signiflcant range of chemical potentials 
and strange quark masses [57]. In the opposite limit where the strange quark 
mass 777.5 is large enough that strange quarks can be ignored, then up and down 
quarks may pair in the 2-flavor superconducting (2SC) phase. Other possible 
condensation patterns are CFL-7^'"'" and CFL-7r°'~, gCFL (gaplcss 

CFL phase), ISC (single-flavor-pairing), CSL (color-spin locked phase), and 
the LOFF (Larkin, Ovchinnikov, Fulde, and Ferrell) crystalline pairing phase, 
depending on the quark flavor densities, the quark chemical potential /iq, and 
electric charge density. For chemical potentials that are of astrophysical inter- 
est, fig < 1000 McV, the color gap is between 50 and 100 MeV, which has a 
signiflcant impact on stellar cooling by neutrino emission. 

2.5 The Surface Photon Luminosity and the Envelope 

The photon luminosity is traditionally expressed as 

= AnR^ ■ as^T^ , (6) 

which defines the effective temperature Te {a^B being the Stefan-Boltzmann 
constant and R the stellar radius). The quantities L, R, and Tc are local 
quantities as measured by an observer at the stellar surface. An external ob- 
server "at inflnity" will measure these quantities red-shifted, i.e., = e^'^L^, 
Too — C^Te, and it!oo = e~'^R, where e^*^ = goo is the time component of the 
metric and is related to the red-shift z by — 1 + z, so that 

= AttRI ■ a,,T^ . (7) 

The luminosity L^, or Loo, is the main output of a cooling calculation, and it 
can equally well be expressed in terms of or T^o- 
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Numerical simulations calculate the time evolution of the internal temperature 
T = T{p,t) and luminosity L = L{p,t) profiles (viewed as functions of the 
density p instead of the radius r) up to an outer boundary pb- This boundary 
is chosen such that at this point the diffusive luminosity I/(pb) is equal to the 
photon luminosity at the surface, i.e., L(pb) = L^, and an envelope model is 
glued as an outer boundary condition. Typically pb is taken as 10^° g cm~^ and 
the envelope is thus a thin layer, of the order of a hundred meters thick, which 
is treated in the plane parallel approximation. Assuming that the thermal 
relaxation time-scale of the envelope is much shorter than the stellar evolution 
time-scale, and that neutrino emission in the envelope is negligible, hydrostatic 
equilibrium and heat transport reduce to ordinary differential equations which, 
with the appropriate physical input, are easily solved. The result is a surface 
temperature Tg = for each given Tb = Ti^p^y). It is usually called a Tb — Tg 
or Tb — Tg relationship. Through Eq. (6), this gives us a relationship between 
L(pb) = and T'(pb) which is the outer boundary condition for the cooling 
code. 

It has been shown in Refs. [58,59] that T^. is actually controlled by a "sensi- 
tivity layer" , where electrons become partially degenerate and ions are in the 
liquid phase. At higher densities the highly degenerate electrons are extremely 
efficient in transporting heat while at lower densities photons take over and 
are also very efficient. The density at which the sensitivity layer is located in- 
creases with increasing T^. This sensitivity layer is hence a throttle and once 
heat has passed through it it can freely flow to the surface and be radiated. 
The layers at densities below the sensitivity layer have no effect at all on the 
thermal evolution of the star, since they are unable to alter the heat flow, 
but the outermost layer, the photosphere, is of course of upmost observational 
importance since it is there that the energy distribution of the emerging flux, 
i.e., the observable spectrum, is determined. 

Gluing an envelope to an interior solution is a standard technique in stellar 
evolution codes. For neutron stars it has two extra advantages: it relieves us 
from solving for hydrostatic equilibrium in the interior, since matter there is 
degenerate, and, most importantly, it allows one to easily include magnetic 
held effects. The magnetic fleld slightly enhances heat transport along it but 
strongly suppresses it in the perpendicular direction, resulting in a highly 
non uniform surface temperature [60]. Assuming that magnetic fleld effects on 
heat transport are negligible at p > pb one keeps spherical symmetry in the 
interior and thus has a unique Tb at pb. For this given uniform Tb one can 
piece together a set of envelope calculations for the various fleld strengths and 
orientations along the stellar surface, corresponding to the assumed magnetic 
fleld structure, and thus obtain a non uniform surface temperature distribution 
Ts{9,4>), in spherical coordinates {9,4>) [61,62]. The effective temperature is 
simply obtained by averaging the locally emerging photon flux F^{9,(j)) = 
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Fig. 6. Two examples of surface temperature distributions induced by the magnetic 
field, in an area preserving projection of the neutron star surface (grey shading, 
shown on the right scale, follows the surface flux instead of the temperature). Left 
panel assumes a dipolar field, with strength 1.2 x 10^^ G at the pole located at 
(6', 4>) = (90°, 90"): for a core temperature of 4.05 x 10^ K it gives Tg = 5.43 x 10^ K 
(see Eq. 8) while the maximum and minimal surface temperatures, at the magnetic 
poles and along the magnetic equator, respectively, arc 

Tmax = 6.70 x 10^ K and 

Tmin = 1-4 X 10^ K. The right panels shows the effect the same dipolar field to which 
a quadrupolar component has been added: this results in Te = 5.31 x 10^ K. This 
particular latter case allows to reproduce the observed ROSAT X-ray pulse profile 
of Geminga (see Figure 6 in [62]) which shows a single very broad pulse while a 
purely dipolar field would result in a double pulse profile (assuming the observer is 
in the direction 9 ~ 90° and emission is isotropic blackbody). Finally, in absence of 
magnetic field, the same internal temperature would result in Tg = 5.54 x 10^ K. 

CSB T^{Q,(t>) over the whole stellar surface 



Two examples of such temperature distributions are illustrated in Figure 6. 
The overall effect on Tc is nevertheless surprisingly small, see, e.g., [61,62] 
and the examples in Figure 6, and a non-magnetic envelope is actually a 
rather good approximation. However, the assumption of spherical symmetry 
at p > pb is questionable and will be discussed in § 7. 

Given that the overall effect of the magnetic field, in the envelope, is not 
very strong, it turns out that the major uncertainty about the envelope is 
its chemical composition. The standard neutron star crust is made of cold 
catalyzed matter, which means ^^Fe at low density (p < 10^ g cm~^). However 
real neutron stars may be dirty and have lighter elements at their surface. 
As was shown in [63] the presence of light elements in the envelope strongly 
enhances heat transport (e.g., the electron thermal conductivity within liquid 
ions of charge Z, in the sensitivity layer, is roughly proportional to 1/2', [64]) 
and results in a significantly higher Te, for the same T^, than in the case 
of a heavy element envelope. Due to pycnonuclear fusion, light elements are 
unlikely to be present at densities above 10^ g cm"^. At very high Tc this 
density is below the sensitivity layer and light elements have little effect, but 
for Te within the observed range (~ 10^ — 10® K) the sensitivity layer is at 
a sufficiently low density so that it can easily be contaminated with light 
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Fig. 7. Relationship between the red-shifted effective temperature and the in- 
terior temperature Tb at the bottom of the envelope assuming various amounts of 
light elements parameterized by 77 = g^i^AMi,/M (AMl is the mass in light ele- 
ments in the envelope, gsu the surface gravity in units of 10^^ cm s~^, and M is the 
star's mass), in the absence of a magnetic field [65]. Also shown are the — Tg°° 
relationships for an envelope of heavy elements with and without the presence of a 
dipolar field of strength of 10^^ G following [66]. Notice that the smaller is AMl 
the lower is the temperature at which its effect is felt. 

elements and the — relationship can be significantly altered. However, if 
only a small amount of light elements is present at the surface their effect will 
only be felt a low T^. We show in Figure 7 the 3], — Tg relationships for various 
amounts of light elements and also, for comparison, the case of a magnetized 
envelope with a 10^^ G dipolar magnetic field. 

2. 6 Comments on minor neutrino emission processes 

The neutrino processes described in § 2.2 are usually the dominant ones but 
in particular cases some minor processes may become important. 

During the first few years of the life of the compact star the surface tem- 
perature is entirely controlled by neutrino emission in its upper layers and 
is independent of what is happening in the core [28,67,68,69]. In this case 
the dominant process is plasmon decay (see, e.g., [48]) which very efficiently 
cools every layer in the the outer crust until the temperature T drops below 
hQp/ks {^p being the electron plasma frequency in the layer) and plasmons 
are exponentially suppressed. This leads to some almost universal temperature 
Te ~ 2.5 X 10^ K (in case of a heavy elements envelope) for all young compact 
stars. Later on the crust will cool by neutrino emission from the electron-ion 
and electron-electron bremsstrahlung processes [48,70] and the PBF process 
from the neutron ^Sq supcrfiuid in the inner crust. This is not very important 
at present time since there are no data about such young compact stars, but 
may become essential in case a compact object is detected within the remnant 
of the supernova SN 1987A. 
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For middle-age stars, in case all processes described in § 2.2 and § 2.3 are 
strongly suppressed by very large gaps then neutrino emission from the crust, 
as described in the previous paragraph, is important. Despite of the small 
amount of mass present in the crust, omission of its neutrino emission may then 
lead to erroneous results, i.e., too warm stars. All numerical results presented 
in this paper of course include these processes. 



3 Some simple analytical solutions 



Assuming simplified physics input one can easily obtain some very illustrative 
analytical solutions to the cooling Eq. (1). Let us write 



^slow _ rj-i8 



(9) 



for the specific heat from degenerate fermions and neutrino emission with only 
"slow" processes or with some of the "fast" processes fisted in Table 1. For 
we write 



usmi 



g reOcr°-^+'* («<!), (10) 



where has been converted into the internal temperature T through an 
envelope model with a power-law dependence. Figure 7 shows that for a heavy 
element envelope or an envelope with a large amount of light elements this 
power-law relation is a good approximation in a wide range of temperatures. 
Typical value for the numerical coefficients in Eqs. (9) and (10) are listed in 
Table 3. 

Due to the much stronger T dependence of compared to L^, at early times 
neutrino emission drives the cooling and when T has sufficiently decreased 
photons will take over. 

(1) During the neutrino cooling era we can neglect in Eq. (1) and, with 
our approximate formulas of Eq. (9), obtain analytical solutions 



Slow u 
coofing 



>: t 



C 
6A^ 



Fast u 
cooling 



y t^-^(^-^], (11) 



47Vf 1^4 



where Tq is the initial temperature at time to = 0. For T <^ Tq, this gives, for 
slow 1/ cooling 



T 



C 
6N' 



1 



and 



(12) 
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Table 3 

Typical numerical coefficients for simplified power-law cooling models 





C 




N' 


S 


a 




erg K"^ 


erg s"^ K"® 


erg 


erg s-i 




high 


1030 


10-32 


10-9 


4 X 10^4 


0.1 


low 


1029 








2 X 10^^ 


0.05 



Comments: 

C: "low" value corresponds to the leptons contribution only, i.e., assuming all baryons are paired with high 
Tc's, while "high" value corresponds to total absence of pairing. 

N^: "low" value corresponds to complete baryon pairing and "high" value to total absence of pairing. 
N^: "high" value corresponds approximately to an inner core of 1 Mq sustaining a direct Urea process with 
nucleons, "low" value of course corresponds to total absence of fast neutrino emission. 
In case of total baryon pairing ~ Aff ~ 0: in this case one has to consider neutrino emission in the 
crystallized crust by the c-ion bremsstrahlung process with a luminosity of the order of LJ"" = AT'^"' . and 
TV" ~ lO-^-* erg s^^ K"*. 

S and a: "high" corresponds to an envelope with a maximum amount of light elements while "low" corre- 
sponds to an envelope made of heavy elements, in both cases with a stellar radius il ~ 12 km. 



and for fast v cooling 



1 



t—^ and T,<xt-s. (13) 



(we have used that a ~ 0). The very small exponent in the Te evolution during 
neutrino cooling is a direct consequence of the strong temperature dependence 
of L,^. The neutrino cooling time scales are also very suggestive: 

T^'°^ = ~ 6 months ■ 



30 



(14) 



and 



_fast 



c 



4 minutes • 



30 



(15) 



and justify the names of "slow" and "fast" neutrino coohng! Notice that 10^ K 
is a typical value for the baryon pairing Tc, and hence, in case of fast neutrino 
cooling, one can expect that a few minutes after the star is born its core may 
become superfluid/superconducting, and the neutrino emission very strongly 
suppressed. 
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(2) During the photon cooling era {Lj S> Li,) one similarly obtains 



^"*^+4^(t^~T^I ' ^^^^ 



where Ti is the temperature at time ti. When t ^ ti and T <^ Ti, we have 

/ C \ ^ 1 ~ 1 

T = — — and Te oc . (17) 

V4q;o / 



Since a <C 1, we see that, during the photon cooUng era, the evolution is 
very sensitive to the nature of the envelope, i.e., a and S, and to changes in 
the specific heat, as induced by pairing. Notice that in case a — one would 
obtain an exponential solution instead of a power law. 

The shift from neutrino to photon cooling: the temperature Tghift at 
which this happens is also easily estimated by equating Lj, to L^, giving for 
slow 1/ cooling 

/ 9 \ 1/6 

rsU^(]^j -lO'K and Te'-IO^K (18) 



and for fast u cooling 

-lO'K and 7;~10^K (19) 



However, for fast u cooling can be significantly reduced by pairing and 
hence T^jft increased. One may also want to obtain an estimate of the age 
4hift at which this happens by using Eq. 12 and 13. It is however now impor- 
tant to keep the a dependence in L^, Eq. 10, and, for small a one obtains, 
for slow u cooling, 11^^^ = C/6S{N^/Sy^/^°'^ and a similar expression for fast 
u cooling. Given that N/S is a very small number, the value of tghift is ex- 
tremely sensitive to the exact value of a, and also of which can be altered 
by several orders of magnitude by pairing. Considering moreover that these 
simple solutions are based on approximate formulas, such results for ighift are 
not very useful in themselves for numerical estimates but show that ighift can 
change significantly by small changes in the interior (through N) and/or sur- 
face (through a) physics. The numerical results presented in the following 
sections show that tguft can vary from about 10*^ yrs up to about 10® years, 
for both slow and fast u coohng. 
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4 The "Minimal Cooling" peiradigni 



The Minimal Cooling paradigm assumes that no enhanced neutrino emission 
is allowed, in particular that no form of "exotic" matter is present in the 
inner core of the star. It is hence a benchmark against which observations 
of cooling neutron stars have to be compared, and a discrepancy between 
the minimal cooling theoretical predictions and data should be considered as 
strong evidence for physics beyond minimal. As such the minimal paradigm 
excludes a prioriihe presence of charged meson condensates, hyperons, and/or 
deconfined quark matter, but it also assumes that the growth of the symmetry 
energy with density is slow enough that the nucleon direct Urea process is 
forbidden. It can be seen as a modern version of the "Standard Cooling" 
scenario but with the important difference that the effects of nucleon pairing 
are wholly taken into account, particularly the strong neutrino emission from 
the formation of Cooper pairs during the pairing phase transition ^ . 

This paradigm is studied in great detail in Ref. [2]. Its tenets immediately 
imply strong constraints on the supernuclear equation of state, resulting in 
stellar radii between 11 to 12 km for a 1.4 Mq mass and between 9 to 10 km 
at the maximum mass (the latter being not so strongly constrained, around 
1.7 up to 2.3 Mq). Neutrino emission from the core comprises the neutron 
and proton branches of the modified Urea processes and the similar, but less 
efficient, n-n, n-p, and p-p bremsstrahlung processes (see Table 1), which are 
inherited from the old "standard" scenario. However, essential and unavoid- 
able within the minimal scenario is the occurrence of nucleon pairing. The 
resulting suppression of both the nucleonic specific heat and the above men- 
tioned neutrino emission processes is well known, but a result is that neutrino 
emission by the PBF process is dominating for most cases of pairing gaps. 
This is illustrated in Figure 8: the left panel clearly shows that for Tc's of the 
order of 10^ K, as in the cases "a" and "T72" , the PBF process is more than 
one order of magnitude more efficient than the unsuppressed modified Urea 
and bremsstrahlung processes when T is below Tg. Once the pairing suppres- 
sion of these last two processes arc taken into account the PBF process can 
dominate by up to two orders of magnitude as shown in the right panel of 
Figure 8. The left panel of Figure 9 compares the effects of various neutron 
— gaps, confirming the results of the previous figure that ~ 10^ K 
provides the strongest PBF neutrino cooling, while the right panel consider 
the effect of the stellar mass, which turns out to be rather weak. Only in case 
of extensive nucleon pairing can we find some small mass dependence. 



^ The motivation for the renaming from "Standard" to "Minimal" is precisely that 
"Standard Cooling" is commonly understood as "Modified Urea Cooling" while, in 
presence of pairing, the dominant process is the PBF process. 
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Fig. 8. Left Panel: neutrino luminosities from the PBF process for four different 
neutron — "'F2 gaps (labeled "a", "b", and "c' as in Figure 4 while the gap 
labeled "T72" is taken from [71]). The lower left panel shows the four gaps' vs. 
the enclosed volume (left scale) and density (right scale) and the upper left panel 
the corresponding PBF luminosities as well as the surface photon luminosity while 
the dashed curve shows, for comparison, the total luminosity from modified Urea 
and nucleon bremsstrahlung processes without pairing suppression. Luminosities are 
calculated assuming an isothermal star so that, e.g., with a temperature of 2 x 10^ 
K, one can see that in case of gap "b" , two regions have unpaired neutrons: the inner 
core at p above ~ 8.5 x 10"*^^ g cm^'^ and a corresponding volume of ~ 1.1 x 10^^ cm'^, 
and the outer core at p below ~ 3 x 10^^ g cm~^ and corresponding volume of about 
~ (11.5 — 9.8) X 10^^ ~ 1.7 X lO-*^^ cm^, while the intermediate region, of volume 
~ 9 X 10^^ cm^ has T slightly below and produces the large Ly by the PBF 
process shown in the upper panel. Right panel: comparison of luminosities from 
various processes during a realistic cooling history: photon ("7")) z^-processes 
("Total z^"), modified Urea and nucleon bremsstrahlung ("MU-Br. z/"), and PBF 
("PBF z^") from n — ^-^2 and p^Sq pairing marked by "n" and "p", respectively. 
The n ^P2 — ^F2 gap is our model "a" which, as shown on the left panel, is amongst 
the most efficient one for the PBF process while the p ^Sq gap is from [72]. All 
results are from Ref. [2]. 

Comparison of the predictions of the Minimal Cooling paradigm with the data 
is shown in Figure 10. Predictions assuming a heavy element or a light element 
envelope are plotted separately and, for each envelope model, the width in the 
predictions is a result of the uncertainty in the size of the nucleon pairing 
gaps. Each wide grey strip encompasses the whole range of results when all 
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Fig. 9. Minimal cooling. Left panel: comparison of cooling trajectories with vanishing 
n — gaps, labeled "0", and the three model gaps "a", "b", and "c" (see 
Figure 4). Right panel: dependence of the cooling trajectory on the stellar mass 
without pairing, where the mass effect is almost indistinguishable, and with only 
proton ^5*0 pairing (dashed curves) and both proton ^Sq and neutron — 
pairing. Results from Ref. [2]. 



gaps shown in Figure 4 are used. The small mass dependence is also included in 
these results. Objects with the best data, shown as boxes in the figure, show a 
remarkably good agreement with the theoretical results of the minimal cooling. 
One may still focus on two stars, PSR 1055-52 and RX J0720.4-3125, which 
may be warmer than predicted and could be cases in which some internal 
heating mechanism is at work. On the other side, PSR 0833-45 ( "Vela" ) and 
PSR 1706-44, may be too cold and require some enhanced neutrino emission. 
However, much stronger cases for the necessity of enhanced neutrino emission 
are the two pulsars PSR J0205+6449 and RX J0007.0+7302 which are clearly 
below any of the predictions of the minimal cooling paradigm. Finally, the 
upper limits on the luminosity of the neutron star which may be present in 
the four SNRs marked as "a", "b", "c", and "d" in this figure doubtlessly 
require enhanced neutrino emission in case any of these SNRs contained a 
neutron star. Similar results has been recently obtained in Ref. [73] which, 
with ad-hoc gaps, showed that even PSR J0205+6449 and RX J0007.0+7302 
could be accommodated within the minimal paradigm. 
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Fig. 10. Comparison of the predictions of the Minimal Model of neutron star cooling 
with the best presently available data. The two dark grey shaded areas correspond 
to models having a heavy elements envelope or an envelope with a maximum amount 
of light elements, as labeled, and the light grey area indicate intermediate trajecto- 
ries corresponding to intermediate amounts of light elements in the envelope. The 
spread of predictions, for each envelope type, corresponds to different assumptions 
about the extend of nucleon (neutron and proton) pairing. Boxes correspond to 
neutron stars where surface thermal emission is clearly detected and which have 
been studied in detail: 1 to 6 are obtained from spectral fits with magnetized hy- 
drogen atmospheres while 7 to 11 are from blackbody fits. These stars are: 1 - RX 
J0822-4247 (in SNR Puppis A), 2 - IE 1207.4-5209 (in SNR PKS 1209-52), 3 - 
PSR 0538+2817, 4 - RX J0002+6246 (in SNR CTB 1), 5 - PSR 1706-44, 6 - PSR 
0833-45 (in SNR "Vela"), 7 - PSR 1055-52, 8 - PSR 06564-14, 9 - PSR 0633-hl748 
("Geminga"), 10 - RX J1856.5-3754, and 11 - RX J0720.4-3125. The next four 
stars, labeled as A, B, C, and D, are barely detected and in case C there is no 
evidence for thermal emission: A - CXO J232327.8+584842 (in SNR Cas A), B - 
PSR J0205+6449 (in SNR 3C58), C - PSR J1124-5916 (in SNR G292.0+1.8), and, 
D - RX J0007.0-F7302 (in SNR CTA 1). The last four data points, a, b, c, and d, are 
from deep observations of four shell SNRs, considered as products of core collapse 
supernovae, in which there is no evidence of any kind for the presence of a compact 
object; they may correspond to very cold neutron stars or isolated black holes: a - ? 
(in SNR G315.4-2.3), b - ? (in SNR G093.3-F6.9), c - ? (in SNR G084.2-0.8), and, 
d - ? (in SNR G127. 1-^0.5). See Ref. [2] for discussion and references. 
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5 Enhanced Cooling 



In the core, almost any chemical composition beyond the one of the minimal 
scenario will open channels for enhanced neutrino emission which possibly 
results in extremely fast cooling of the neutron star. The left panel of Figure 
11 illustrates the effect of the direct Urea process from nucleons in a model 
where this process is allowed for masses above 1.35 Mq (this specific critical 
mass is of course very model dependent). Notice that the 1.4 Mq star has 
an inner "pit" where the direct Urea process is occurring with a mass of only 
0.038 Mq , while the 1.7 Mq star's "pit" is above 1 Mq . After the early phase, 
at ages a few tens of years, the rapid temperature drop in fast cooling stars 
is due to the finite thermal relaxation time of the crust and heavier stars, 
having thinner crusts, relax faster. The right panel of this figure compares 
the impact of the direct Urea process on stellar cooling with other enhanced 
cooling processes, which originate from pion condensation, kaon condensation 
and the direct quark Urea process. Differences between these models come as 
much from the various efficiencies for neutrino emission as from the various 
critical densities at which enhanced emission becomes allowed. Comparison 
with the data plotted here shows that any of these fast coohng processes 
results in stars with temperature clearly incompatible with observations. 

This naive picture is strongly affected when baryon pairing is taken into ac- 
count [30,74,75]. As discussed in § 2.3 the development of an energy gap in- 
duces a suppression of the neutrino emission for any process in which the 
paired component participates. As can be seen from the left panel in Fig- 
ure 12, neutron superfluidity (^P2 — '^^2 in the present example with a small 
superfluid gap of 0.3 MeV [5]) delays cooling signiflcantly and moves the the- 
oretical temperatures right up to the region where most observed data are 
concentrated. Similar results are obtained with all enhanced cooling scenarios 
as, e.g., nucleon direct Urea [75], vr" condensates [76,77], hyperons [78,79,80] 
or hyperons with deconflned quarks [80] . With a high enough value of Tc the 
enhanced neutrino emission may not even have time to act and result in stellar 
temperatures as high as the ones obtained within the minimal scenario [75]. 

The right panel of Figure 12 illustrates several of these considerations [80]. 
This panel shows cooling curves for models with only nucleons (thin contin- 
uous curves) and with nucleons and an hybrid phase of nuclear+deconflned 
quark matter (dashed curves and dotted curves) and various hypotheses about 
pairing. Nucleon and quark direct Urea processes are allowed, except in the 
case marked as "Mc" where they are arbitrarily turned off for illustration. Four 
different neutron — ^F2 gaps are considered, labeled as "z" (zero gap), and 
"a" , "b" , and "c" according to Figure 4. For models including quark matter in 
the inner core, five quark gaps are also considered, labeled as "Z" (zero gap), 
"A", "B","C", and "D" corresponding, to maximum gap sizes of 0.1, 1.0, 10, 
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Fig. 11. Left panel: influence of the direct Urea process on the cooHng of stars of 
various masses. (Figure from Ref. [75].) Right panel: influence of different enhanced 
neutrino emission processes on the cooling of neutron stars of mass M = 1.4 M©, 
except for the kaon-condensed model whose mass is M = I.SMq. (The 1.4 
model is not dense enough to support a kaon condensation). (Figure from Ref. [32].) 
Upper limit on PSR 1929+10 thermal luminosity is from Ref. [82]; see Figure 10 for 
identification of other data. 

and 100 MeV, respectively (all flavors and color of quarks are assumed to 
have the same gap for simplicity). Considering models without quarks matter 
(continuous lines), one sees that models with a large neutron — gap 
"c" have almost the same evolution as the model with no direct Urea process, 
"Mc", and are all compatible with the highest measured neutron stars tem- 
peratures. Models with neutron gap "b" could explain, as in the left panel, the 
intermediate temperature neutron stars, in which case the warmer one should 
be understood as having lower masses and no direct Urea process allowed 
(model "Mc"). Models with smaller neutron gaps produce too low tempera- 
tures but would correspond to the four upper limits (data points labeled as 
"a", "b", "c", and "d" in Fig 10). Finally, hybrid models with quarks (dashed 
and dotted lines) give almost indistinguishable results when the quark gaps 
are large enough (curves labeled "A" , "B" , "C" , and "D" ) and it is only in 
case of a vanishingly small quark gap that their presence is noticeable. 

As can bee seen from this brief presentation, we have an embarrassingly large 
number of possible scenarios in case enhanced neutrino cooling is occurring. 
Most models have to include in an essential way the various possible pair- 
ing gaps, about which very little is known at high densities, to reconcile en- 
hanced neutrino cooling with data. Nevertheless there are possibilities to have 
a smooth transition from slow to fast neutrino cooling through strong medium 
effects and without invoking pairing [22,81], or with a significant amount of 
internal heating as discussed in § 8. Obviously, much work is still needed to 
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Fig. 12. Left panel: influence of superfluid — neutrons on the cooling curves 
labeled GJoq and G300 (with kaon condensation) in Figure 11 for various neutron 
star masses. (Figure from Ref. [32].) Right panel: effect of the neutron and quark 
gap size. See text for description and discussion. (Figure from Ref. [80].) The large 
differences in luminosities at late times between the two panels are due to small 
differences in treatment of the envelope and specific heat, as discussed after Eq. 17. 
See Figure 10 and 11 for identiflcation of data. 



determine which scenarios are the more plausible. 



6 Cooling of Strange Queirk Steirs 



6.1 Strange Stars with Nuclear Crusts 



If strange quark matter were in fact the true ground state of the strong inter- 
action (see § 2.1), new classes of compact stars should exist which range from 
dense strange stars to strange dwarf stars [83,84]. They would form distinct 
and disconnected branches of compact stars and are not part of the continuum 
of equilibrium configurations that include ordinary white dwarfs and neutron 
stars [85]. Figure 13 shows the cooling behavior of such strange stars with 
and without nuclear crusts. One sees that not even the thickest, theoretically 
possible crust (inner crust density equal to neutron drip density) does prevent 
strange stars from cooling very rapidly [86]. If one treats the quarks as su- 
perfluid particles, assuming a small density-independent gap of just 0.1 MeV, 
radiation of neutrinos is greatly reduced and the stellar cooling behavior is 
rather similar to that of conventional neutron stars. The different condensa- 
tion patterns of color superconducting quark matter, discussed in § 2.4, do not 
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Fig. 13. Cooling behavior of strange stars with maximaUy thick crusts and with 
very thin crusts, i.e., with only an envelope (Figure from Ref. [32].) See Figure 10 
and 11 for identification of data. 



alter this picture dramatically. In the CFL phase, for instance, all quarks have 
very large gaps, A ^ T, so that both and are so strongly reduced that 
quark matter that may exist in the center of a neutron star would be rendered 
invisible. Notice that, in spite of the CFL phase being charge neutral by itself 
in bulk [87] , a strange star in the CFL can still support a nuclear crust when 
surface boundary conditions are properly taken into account [88]. The cool- 
ing behavior of such stars would then be determined by the nuclear matter 
surrounding the CFL quark matter core. This conclusion is strengthened by 
the studies performed in Refs. [89,90]. The cooling behavior of compact stars 
with 2SC quark matter in their cores is simplified by the fact that up and 
down quarks may pair with a gap A ~ 100 MeV which is orders of magnitude 
larger than the stellar temperature, ^ 1 MeV, and are therefore inert with 
respect to the star's temperature evolution. For 2SC quark matter, however, 
there exist also quark pairing channels that lead to weak pairing with gaps 
on the order of several keV to about 1 MeV, which is on the same order of 
magnitude as the star's temperature. These quarks may thus not pair but, in- 
stead, radiate neutrinos rapidly via the quark direct Urea process (Table 2). If 
this is the case, the 2SC quark matter core would cool rapidly and determine 
the cooling history of the star [55,91]. Naturally the cooling behavior of such 
stars depends rather sensitively on the value of the superfluid gap [91]. 
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6.2 Strange Stars with a Bare Quark Surface 

Even though it is in principle possible to cover the strange quark matter by 
a layer of normal nuclear matter, as assumed in the previous subsection, it is 
not clear at all how STich a nuclear crust may be formed. When a strange star 
is born its internal temperature is likely of the order of a few times 10^^ K, 
as in a standard proto-neutron star formed in core-collapse, and the resulting 
neutrino flux is so high that it should easily be able to expel all baryonic 
matter surrounding the quark matter [92] . As a result a new-born strange star 
is most certainly bare, i.e., its surface consists directly of quark matter with 
no nuclear component above it. Such a strange star remains bare as long as 
its temperature is above ~ 3 x 10^ K [93]. 

If this is the case, the thermal evolution of a bare strange star would be rad- 
ically different from the one of a neutron star or a strange star with a crust. 
The surface density of the quark matter is about 4 — 8 x 10^^ g cm~^ and the 
resulting plasma frequency of the order of 20 MeV [18] which may lead one to 
believe that a bare strange star is unable to emit thermal radiation once it's 
temperature dropped below about 10^° K, i.e., a few seconds after its birth: 
a bare strange star would be a silver sphere instead of a blackbody emitter. 
However this naive picture neglects the presence of electrons within the quark 
matter which are bound to the quark matter by Coulomb forces. As such, the 
electrons slightly leak out of the quark matter, by a distance which is roughly 
given by their Debye screening length, producing thus an electrosphere with 
a thickness of a few thousand Fermis [94] . As noted by Usov [95] the result- 
ing electrostatic field, of the order of 5 x 10^'' V cm~^ [18], is well above the 
Schwinger critical field of QED and can induce copious production of e~ — 
pairs at the surface. However, this pair emission does not occur in vacuum, 
but within the electrosphere and it is ultimately limited by electron degener- 
acy. Nevertheless the resulting luminosity is enormous as long as the surface 
temperature is above ~ 10^ K [96]. Pairs outflowing from the stellar surface 
mostly annihilate into photons in the vicinity of the strange star [96] resulting 
in a fireball of e^, and 7's. At luminosities > 10'^^ erg s^^ this plasma is 
optically thick and produces a blackbody spectrum while at decreasing lu- 
minosities the spectrum progressively evolves into a very wide annihilation 
line [97]. At lower temperatures, when pair production become suppressed by 
electron degeneracy, photons can still be emitted by — e+ bremsstrahlung 
within the electrosphere [98] and result in significant luminosities. 

The thermal evolution of a bare strange star with the Usov- Schwinger pair 

production mechanism has been studied numerically in [99] and results are 
shown in Figure 14. Uncertainties about the pairing state of quark matter are 
taken into account schematically in the three scenarios marked in the figure 
as "A", "B", and "C", in which the gaps suppress neutrino emission with 
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Fig. 14. Cooling behavior of a bare strange star, showing the thermal luminosity 
from the Usov-Schwingcr pair-production mechanism vs. age. In the left panel it is 
assumed that heat, within the star, flows to the surface by diffusion while in the 
right panel convection is allowed resulting in much stronger heat flow and higher 
thermal luminosities. In the dotted curves no quark superconductivity is included 
while the continuous curves labeled A, B, &; C have quark color superconductivity 
taken into account for three possible scenarios (see text for details). (Prom Ref. 



increasing efficiency, leading thus to slower cooling of the star. These results 
show that bare strange stars can produce extremely high thermal, luminosi- 
ties, well above the Eddington limit, over extended periods of time, but with 
considerable uncertainty due to the precise phase of color superconductivity 
present at these "low" astrophysical densities. Notice, however, that the pho- 
ton emission from e" — bremsstrahlung in the electrosphere [98] was not 
included in these models, and will significantly change the evolution when the 
thermal luminosity has dropped below 10^° erg s~^. 



7 Magnetic field effects in the crust 

All the heat stored in the core of the neutron star and eventually irradiated 
away from its surface by photons has to be transported through the crust. 
In the absence of rotation and magnetic field this transport in the crust is 
spherically symmetric. While the effects of rotation are quite small even for 
millisecond pulsars, the presence of magnetic fields may cause significant devi- 
ations from the spherical symmetry of the transport processes, even for quite 
"standard" field strength of ~ lO^^G. Due to the classical Larmor rotation of 
electrons, a magnetic field causes anisotropy of the heat fiux since the heat 
conductivity becomes a tensor whose components perpendicular, k^, and par- 
allel, K\\, to the field lines become 



[99].) 




and K\\ — Ko 



(20) 



29 



Fig. 15. Temperature distribution and magnetic field lines in the crust for a "core" 
(left panel) and a "crustal" (central panel) dipolar field (from Ref.[100]), and sur- 
face temperature distribution resulting from the "crustal" field (right panel). The 
thickness of the crust in left and central panels has been stretched by a factor 5 for 
clarity and the two bars show the temperature scales in units of Tcore = 

10^ K. The 

"core" field gives a surface temperature distribution very similar to the left panel 
of Figure 6. For the "core" filed on has Tg = 1.15 x lO^K while for the "crustal" 
field Te = 0.89 x lO^K. In both cases, the field strength is = 3 x lO^^ q (xhg 
neutron star considered here has a 1.4 Mq mass and a radius of 11 km.) 

where is the conductivity in absence of magnetic field, Qb the cyclotron 
frequency of the electrons and r their coUisional relaxation time. In the neutron 
star crust Qbt ^ 1 is easily realized and, as a result, heat flows preferentially 
along the magnetic field lines [100]. This effect is moreover amplified at the 
surface by the well known non-isotropy in the envelope as described in § 2.5. 

The non-isothermality of the subjacent crust depends strongly on the internal 
geometry of the field [100]. Assuming a dipolar field structure, outside the star 
the radial dependence of the field is uniquely determined, oc r~^, while inside 
the star it depends on the location of the electric currents, through Ampere's 
law. Assuming currents are exclusively located in the core, "core field", the 
dependence also applies to the crust while if they are exclusively located in the 
crust and the field does not penetrates the core, "crustal field", the field topol- 
ogy in the crust is radically distinct: for the same external field in the latter 
case field lines in the crust are squeezed into the crustal shell and have hence 
a very large meridional component. This meridional component of the crustal 
field inhibits the radial flow of heat which is hence redirected preferentially 
toward the magnetic poles. This difference is illustrated in Figure 15. 

The drastic difference in the crustal temperature distribution for the different 
field structures which are characterized by the same dipolar field structure and 
strength outside the neutron star, causes significant differences in the surface 
temperature distribution which will have several observational consequences: 

1. A non- uniform surface temperature induces a modulation of the observed 
soft X-ray thermal emission [61]. The stronger channeling of the heat flow 
toward the polar regions in the case of a crustal field will result in larger 
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amplitude in the pulse profile. 



2. The differences in the photon luminosities for a core or a crustal field will also 
affect the long term cooling of neutron stars. A neutron star having a magnetic 
field confined to its crust will stay warmer for a longer time, due to its lower 
photon luminosity, compared to a neutron star with a field penetrating its 
core. 

2. This may open a way to study the internal geometry of the magnetic field. 
We may be able to distinguish between field geometries where the currents 
are essentially localized in the crust or in the core, or even detect the imprint 
of the presence of a toroidal field [101] 



8 Heating Mechanisms 



The expression "heating mechanism" refers to the term in Eq. 1 and 

generically encompasses all possible dissipative processes which will inject heat 
into the star by tapping into various forms of energy: magnetic (§ 8.1), rota- 
tional or chemical (§8.2). 



8.1 Magnetic Field Decay and Joule Heating 



Given that most neutron stars have strong magnetic fields, magnetic energy 
is a natural reservoir from which to extract heat by the Joule effect from 
the decaying electric currents. Assuming a uniform internal field of strength 
B — 10^^ Bis G, one can roughly estimate an amount 

£^mag ~ 1^ X ^TtR^ - 2 X 10^^ ^.g (21) 



of stored magnetic energy. With a field decay time scale r = lO^Te yrs this 
gives us an equivalent "magnetic heating luminosity" 

^mag ^ ^ 6 X 10^^ ^ erg s'^ (22) 



This simple estimate indicates that decay of a standard magnetic field can 
alter the neutron star thermal evolution at late (> 10^ yr) times. In the 
case of a magnetar, ultra-strong fields and potential channels for fast decay 
(i.e., r <C 10^ yrs) can nevertheless lead to significant heating in young stars 
[102,103,104]. 
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We will now describe some numerical results for the simplest case of a stan- 
dard magnetic field, following the Ref. [105] which considered the case where 
magnetic energy is converted into heat through Ohmic decay of the currents 
supporting the field. This Joule heating contributes to the source term H in 
Eq. (1), where the heat production by field decay per unit of (proper) time 
and volume is given by j^/o", j being the electric current density and a the 
scalar electric conductivity. Joule heating is specially effective in the crust, 
which never becomes superconductive. Therefore, Ref. [105] considered only 
magnetic fields which are, along with the currents supporting them, confined 
to the neutron star crust and some results are shown in Figure 16 (see also 
Ref. [106].) The amount of heat locally released by Joule heating is determined 
by the strength of the field and its decay rate at that time. Most of the mag- 
netic energy is dissipated at early times when it has almost no effect because 
the star's thermal energy is still too large. It is only when the star has be- 
come sufficiently cold that I/mag becomes significant, unless the field strength 
is of magnetar size. As illustrated in Figure 16, scenarios with a weak field 
supported by currents located at low density, i.e., decaying fast while the star 
is still hot, lead to virtually no observable effect. On the contrary, scenarios 
with strong field supported by currents at high densities keep a large amount 
of magnetic energy to be dissipated on long time scales, the best cases being 
fast cooling stars which result in the highest temperatures at ages ~ 10^ — 10^° 
years. 

The above description is based on a linear evolution of the field but non-linear 
effects caused by the Hall-drift in the highly magnetized crust may lead to an 
instability, which transfers magnetic energy rapidly from large scale structures 
into much smaller ones [107,108]. In that way, close to the NS surface small 
scale (~ 0.1 — 1 km) field structures, as necessary for the onset of the pulsar 
mechanism, will be created which simultaneously are sources of accelerated 
(^Haii ~ lO^yrs) Joule heating [109]. Detailed numerical modehng of the effect 
of this instability on the cooling of the pulsar remains to be performed. 



8.2 Dissipative motion of vortex lattices and readjustment to equilibrium 

The dissipative motion of vortex lattices and the rotational readjustment of 
the stellar equilibrium structure are other key processes that contribute to the 
heating of neutron stars. The relative effectiveness of these heating processes, 
however, varies significantly from one process to another, and depends sensi- 
tively on the value of the rotational parameter K, which enters the power-law 
relation Q{t) = —K{t)Q"'{t) characterizing the spin-down of pulsars. Figure 
17 shows the heating rate H (which enters in Eq. (1)) computed for several 
competing, internal heating processes and parameter sets. These are [110]: A) 
vortex pinning at nuclei for the Epstein-Baym parameter set (EB- pinning) ; 
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Fig. 16. Thermal (upper panels) and magnetic (lower panels) evolutions of neutron 
star undergoing slow (left panels) and fast (right panels) neutrino cooling, with Joule 
heating from magnetic field confined to the stellar crust. The three representative 
case of initial field strengths cover a wide range of model parameters: Bq = 10^^, 
10^^, and 10^^ G have their supporting currents initially concentrated at densities 
of po ~ 10^^, 10^^, and 10^^ gm cm~^, respectively. (Adapted from Ref. [105].) 

B) vortex pinning at nuclei for the parameter set of Pizzochero et al.; C) co- 
rotating crust; D) co-rotating core; E) crust cracking; and F) chemical heating. 
Chemical heating occurs because of the changing chemical composition with 
stellar frequency. The matter would maintain chemical equilibrium if the re- 
laxation timescales for the weak reaction processes were small compared to 
the timescale of rotational evolution. These timescales, however, were found 
to be comparable and the stellar composition therefore departs from chemical 
equilibrium, which modifies the reaction rates [111] and leads generally to a 
net conversion of chemical energy into thermal energy [112,113]. In general, 
internal heating leads to enhanced stellar surface temperatures, which can be 
very pronounced for rotating neutron stars of middle {K ~ 10~^^ s) and old 
{K ~ 10^^^ s) ages [110] or millisecond pulsars [114]. As shown in [110,115] 
(see also [116]) , the increase in surface temperature caused by the internal 
heating due to thermal creep of pinned vortices and the outward motion of 
proton vortices in neutron stars cooling rapidly via one of the enhanced mech- 
anisms discussed in § 5 would lead to a good agreement with the observed 
data. 
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Fig. 17. Heating rate H as a function of internal star temperature for the heating 
scenarios (A)-(F) hsted in the text. The left panel corresponds to ^ lO"-*^^ s, 
the right panel to if ~ 10"^^ s. The photon and neutrino luminosities of standard 
cooling and enhanced cooling models (L-y, L^,^ gtd., L^,^ enh.) are shown for comparison. 
(Figure from Ref. [110].) 

9 Cooling Neutron Stairs in Soft X-Ray Transients 

Neutron star undergoing accretion from a companion in a close binary system 
can also give us information about the state of dense matter at supranuclear 
densities. In case of continuous accretion, the bulk of the observable X-ray 
emission is dominated by the release of gravitational energy when matter hits 
the neutron star surface and by frictional heating in the inner part of the 
accretion disc, in case where such a disc exists. Nuclear energy is released 
during X-ray bursts, and superbursts, which may teach us about the thermal 
response of the stellar core [117]. 

A very interesting class of accreting neutron stars arc the one in the so-called 
Soft X-Ray Transients ( "SXRTs" ) which undergo recurrent surges of activity 
separated by long phases of relative "quiescence" during which accretion very 
likely does not occur. The substantial X-ray brightening observed in outburst 
(Lq ~ 10'^''"^^ erg s^^) together with type I bursts unambiguously indicate 
that episodes of intense accretion occur onto a stellar surface, i.e., a neutron 
star, and not into a black hole. In quiescence the faint X-ray emission has a 
luminosity in the order of Lq ~ iQ^^-sa g-i^ j comparable to a standard 
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isolated cooling neutron star. Typically, the duration of an outburst, to, is 
considerably shorter than the recurrence time between two outbursts, t,.. In 
principle, the effects of accretion are the same as in case it is continuous but 
during the quiescence phase we are most certainly detecting thermal emission 
from the cooling of the neutron star heated up by the accretion. Given the 
high internal temperature of the star, the heat released in the upper layers 
from the accretion and the bursts flows back to the surface and is radiated 
away because of the large temperature gradient in the envelope. Nevertheless, 
the star stays hot due to energy release from non-equilibrium process in the 
interior: this heat can be stored in the stellar interior and slowly released when 
accretion stops. 

Non equilibrium processes certainly occur in the crust of an accreting neu- 
tron star: thermonuclear processes at the surface will produce iron-peak nu- 
clei which are then pushed to higher densities. In their way, these nuclei will 

undergo electron capture, neutron emission and absorption and eventually 
pycnonuclear fusions [118]. Overall, about 1.5 MeV is released as heat for 
each accreted baryon. This energy is enough to explain the observed quiescent 
(thermal) luminosities Lq in most cases [119]: 

L^-fx Qnuc^^ ~ / X 6 X io32_^£^^ ergs'S (23) 

^ 1.5MeV IQ-iiM© yr"! ^ ' ^ ^ 



where < M > is the time average of the accretion rate, Qnuc the energy released 
deep in the crust per accreted baryon, and / the fraction of Qnuc which is stored 
in the stellar interior, i.e., not lost by neutrino emission. The luminosity 
during an accretion outburst can be estimated as Lq ~ {AM/to){GM/R) 
where AM is the mass accreted during the outburst. Writing < M>~ AM/t,- 
one then obtains, following [119], 

Loto GM m„ 200 

~ ~ . (24) 

Lqtj- R fQnnc f 



Notice that Lq/ Lq is independent of the source's distance, which is often poorly 
constrained, while to and tr can in principle be directly obtained by monitoring 
the source for a long enough time. 

Detailed numerical calculations [120,121] confirmed this simple and elegant 
picture and are compared with data in Figure 18 (see also [122,123]). Results 
in this figure are based on a dense matter model in which the direct Urea 
process is allowed at high densities and controlled by neutron pairing. The 
two stars Aql X-1 and 4U 1608-52 show very good agreement with the sim- 
ple model and the numerical results with / — 1, i.e., with negligible neutrino 
losses and thus very probably are low mass neutron stars. The star in Cen X-4 
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Fig. 18. Quiescent to outburst luminosity ratio plotted versus the recurrence time 
over the outburst time. The filled squares represent observed values. Estimates 
for the efficiency of heat storage / = 1 and 0.01, Eq. (23) (but with minor GR 
corrections), arc shown for a 1.4 Mq star (dashed-dot) and 1.8 Mq star (dashed). 
Dots connected by lines show numerical results (see text ) for different neutron star 
masses as labeled: the 1.4 Mq case has no enhanced neutrino emission, i.e., / ^ 1, 
while the other two cases at 1.7 and 1.8 Mq have increasingly more efficient neutrino 
emission, i.e., / ^ 1 (the precise value of the neutron star masses are, however, very 
model dependent). Adapted from [120], with new data for SAX J1808. 4-3658 from 
[124]. 

is noticeably off and seems to require very strong neutrino losses. However, 
it must be emphasized that in 35 years since its discovery it has undergone 
only two bursts (in 1969 and 1979) and hence its recurrence time is highly 
uncertain: it is not impossible that this point should be plotted much further 
to the right, as indicated by the arrow in the figure. The most interesting and 
intriguing object is certainly SAX J1808. 4-3658, dubbed as "The accreting 
pulsar", which provided the first strong observational evidence for the evolu- 
tionary scenario that millisecond pulsars owe their short rotational period to 
a long phase of accretion in a low mass X-ray binary system. The very low 
upper limit on its thermal quiescent luminosity Lq [124] requires, within the 
present scenario, very strong neutrino emission (even beyond what has been 
considered in [120] and shown in Figure 18, i.e., an inner core where enhanced 
neutrino emission is not suppressed at all by neutron pairing) and hence a 
large mass which is in perfect agreement with the idea that recycling of a neu- 
tron star to a millisecond rotation period by accretion requires a mass transfer 
of about 0.2 to 0.5 Mq [125]. Notice that the first precise measurement of the 
mass of a millisecond pulsar, PSR J 1909-3744, was recently announced, giv- 
ing a value of 1.438 ± 0.024 Mq [126], at the upper edge of the mass range 
of binary pulsars [127]. (The latter, which are not millisecond pulsars, come 
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from short-lived double massive star systems and did not have time to ac- 
creted a significant amount of mass.) The measured mass of PSR J1909-3744 
demonstrates that spinning-up of a neutron star to millisecond period does 
not requires accretion of much more that 0.2 Mq . Other mass estimates of bi- 
nary millisecond pulsars [128], though not as accurate, do indicate that much 
larger masses are reachable, e.g., PSR J0751-I-1807, which has a mass lower 
limit of 1.6 M0 . 

Another, and still poorly explored, aspect of these objects is their short time- 
scale thermal response to the accretion phases. This is an aspect of the problem 
from which extremely important information about both the structure of the 
neutron star crust and the thermal state of their core can be obtained [121,129], 
and about which intriguing observational results are being found (see, e.g., 
[130]). 

10 Epilogue and Conclusions 

Confrontation of the predictions of the Minimal Model with data in § 4 showed 
a reasonably good agreement with presently available data on isolated cool- 
ing neutron star, i.e., no compelling evidence for the occurrence of enhanced 
neutrino emission is found. Nevertheless, most models of dense matter predict 
the presence of some form of exotic matter at high enough density, or at least 
a high proton fraction, which allow enhanced neutrino emission. Are these 
"non-minimal" scenarios wrong or are we missing cooling neutron stars ? 

Before hasting into conclusions, one must first emphasize that the fact the 
results of the Minimal Cooling scenario arc compatible with the data does not 
mean that the other scenarios are ruled out (unless one is a fervent adept of the 
"Occam's Razor Principle"). A scientific paradigm is refuted when it is incom- 
patible with experimental facts, not because its competitors are compatible. 
Alternative scenarios, "exotic" ones or simply with a high proton fractions, 
whose chemical composition allows enhanced neutrino emission do not nec- 
essarily imply fast cooling since the neutrino emission can be suppressed by 
pairing we showed in § 5. With large enough gaps these scenarios can be com- 
patible with the data as illustrated in Figure 12. Moreover, enhanced neutrino 
emission is expected to be allowed, if at all, only at high enough densities, 
and hence for sufficiently massive neutron stars, the critical neutron star mass 
depending on the particular model. 

We have no direct information at all about the masses of the observed nearby 
cooling neutron stars. Many measured masses of radio pulsars in binary sys- 
tems are between 1.25 to 1.44 Mq [127,131]. Exceptions to this are the mass 
of the binary neutron star J0751+1807 which is in the range 2.ltg;5 Mq [128], 
and the companion of PSR J1756-2251, with a mass of 1.18l[j:[]^ Mq [132]. 
It is not clear if the mass range of 1.25 to 1.44 Mq also applies to isolated 
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neutron stars. Evolutionary models of massive stars [133] predict a bimodal 
distribution for the neutron star masses, with the two peaks being between 
1.2 — 1.3 Mq for progenitor main sequences masses below 19 Mq (which will 
burn carbon in a convective core) and between 1.7 — 1.8 Mq for progenitor 
main sequence masses above 19 Mq (where carbon burning will occur in a 
radiative core). The actual value of 19 Af© for the main sequence bifurcation 
mass may not be very accurate, but the double peaked neutron star mass 
distribution is certainly real. Given the initial mass function for heavy stars, 
one could thus expects that massive (e.g., > 1.5 Mq ) isolated neutron stars 
be very rare. 

A neutron star whose temperature is in disagreement with the Minimal Cool- 
ing predictions must be of course cold, but moreover young and very likely 
massive, hence probably uncommon. Therefore, such a star can be expected 
to be found far away from the Sun and its detection will be very difficult since 
it's weak thermal spectrum will be strongly absorbed. Notice that all young 
(< a few times 10^ yrs) neutron stars plotted in Figure 10 are at a distance 
superior to 1 kpc except for the Vela pulsar. In particular, the two most con- 
spicuous young and cold objects, PSR J0205+6449 (in SNR 3C58) and RX 
J0007.0+7302 (in SNR CTA 1) are at distances of 2.6 ± 0.3 kpc and 1.4 ± 0.3 
kpc, respectively. Moreover, as emphasized recently in [134], about 80% of the 
nearby young neutron stars come from the Gould Belt: this local structure of 
young stellar associations forms a ring of about 1 kpc surrounding the Sun 
and has an estimated age of 30-50 Myrs. If the Gould Belt's progenitors of the 
nearby young cooling neutron stars were formed when the belt was formed, 
which is likely for many of them, and core collapsed about one Myrs ago, they 
were massive stars not much above 10 Mq since heavier stars died well be- 
fore. From the results of [133] one can safely conclude that most nearby young 
cooling neutrons stars have masses not exceeding 1.4 Mq-. it is not surprising 
that none of them shows clear evidence of enhanced neutrino cooling. 

The smoking guns in favor of enhanced neutrino cooling are certainly, to date, 
the two isolated pulsars PSR J0205+6449 and RX J0007.0+7302 (see § 4) and 
the accreting pulsar SAX J1808. 4-3658 (see § 9). Since evolutionary scenarios 
for recycling pulsar to millisecond periods require accretion of about 0.2 to 
0.5 Mq [125], SAX J1808. 4-3658 is a natural candidate for a neutron star 
with an inner core beyond minimal. Possibly more such objects will be found 
in the SNRs studied recently in [135] (whose present luminosity upper limit 
are marked as "a", "b", "c", and "d" in Figure 10) in case a neutron star 
is detected and will provide us with definite evidence for the occurrence of 
enhanced neutrino emission. Finally, one must entertain the idea that core 
collapse supernovae fail to produce heavy neutron stars, sending them into 
a black hole despite of them being allowed by the equation of state, and 
that heavy neutron stars can only be found in binary systems, such as SAX 
J1808.4-3658, after substantial mass accretion occurred. 
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